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POSSIBILITY  OF  DISSOCIATION  REACTIONS 
OCCURRING  IN  POLYATOMIC  MOLECULES 
DURING  RAPID  COOLING  IN  AN  EXPANDING 
SUPERSONIC  FLOW 


P.  S.  Zaslonko,  S.  M.  Kcgarko,  and 
Yu.  V.  Chirikov 

(Moscow) 


Studied  is  the  possibility  of  the 
occurrence  of  the  dissociation  reaction 
of  N20  and  the  exchange  reaction  of  No0+ 

+C0-»-N2+C02  in  rapid  cooling  of  a  gas  mixture 

as  a  result  of  supersonic  expansion  in  a 
nozzle.  Balance  equation  for  vibrational 
energy  and  the  kinetics  of  the  chemical  re¬ 
action  under  nonequilibrium  conditions  are 
solved  numerically.  The  values  of  param¬ 
eters  necessary  for  experimental  observations 
of  the  studied  phenomenon  are  found. 


to  the 

existing 

repre- 


The  rapid  exchange  of  vibrational  quantum  usually  leads 

establishment  of  quasi-equilibrium  vibrational  distributions 

W  VT  W  VT 

in  such  a  course  of  time  t  that  t  <t<x  ,  where  t  and  r 
sent  the  times  of  V-V  and  V-T  processes,  respectively 


For  example,  in  a  system  of  oscillators  A  and  B  when  V-V 
equilibrium  is  reached,  the  following  relationship  between 
temperatures  and  Tg  is  established  [1,  2] 
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efl/rB  =  oa*ta  -  (qb  -  e^/r 


(l) 


where  8A  and  6g  are  the  characteristic  temperatures  of  oscillators 
A  and  B,  respectively,  T  -  gas  temperature. 

Under  certain  conditions  (0g>0A,  T^,  Tg>>T)  there  may  occur 
a  significant  vibration  pumping  of  oscillators  with  a  lower  quantum 
energy  from  oscillators  with  a  higher  quantum  energy.  In  a  system 
of  anharmonic  oscillators  under  nonequilibrium  conditions  the  pump¬ 
ing  of  energy  to  higher  levels  may  lead  to  a  significant  increase 
in  the  ‘'temperature"  of  these  levels.  The  population  distribution 
with  respect  to  the  ievels  of  an  anharmonic  oscillator  (generally 
the  term  "Trinor  distribution"  is  used)  is  described  by  the 
equation  [3,  4] 


AT„  =  ~  exp  i)  ^£]}  =h  -J-j. 


(2) 


Here  Nn  is  the  concentration  of  particles  at  level  n,  Q  -  the 
statistical  sum  of  the  anharmonic  oscillator  and  its  anharmonicity 
(°K),  respectively,  T^  -  effective  vibration  temperature  for 
levels  1-0.  Since  the  rates  of  many  chemical  reactions  are  determined 
by  the  population  of  rather  high  vibrational  levels,  it  is 
interesting  to  examine  the  effect  of  energy  pumping  under  non¬ 
equilibrium  conditions  on  the  occurrence  of  certain  model  reactions. 
The  dissociation  and  exchange  reactions  are  convenient  in  this 
case.  The  constant  of  the  dissociation  rate  is  determined  by  the 
degree  of  vibrational  excitation  of  the  dissociating  molecules.  The 
dependence  between  the  rate  of  the  exchange  reactions  and  vibrational 
excitation  of  the  original  reagent  has  apparent^  not  been 
sufficiently  studied,  although  the  series  of  studies  of  Bauer  and 
co-authors  (see,  for  example,  [5])  definitely  showed  that  vibrational 
excitation  played  a  deciding  role  in  a  great  number  of  exchange 
reactions.  The  experiments  in  HI  and  Di  molecular  beams  showed 
L 6 ]  that  these  molecules,  even  at  a  very  high  translational  energy 
(from  20  to  109  Cal/mole  in  a  mass  center  system,  i.e.,  high 
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activation  energies)  do  not  enter  into  the  exchange  reaction  where 
HD  is  formed  unless  they  possess  sufficient  vibrational  energy. 


Thus,  if  we  are  limited  to  merely  the  gas  dynamic  methods  of 
cooling  (for  example,  the  flow  in  a  nozz?e),  then  the  factors  which 
cause  the  reactions  to  occur  at  low  temperatures  are  as  follows : 

1)  nonrescnance  vibrational  pumping  of  "working"  molecules 
from  inert  biatomic  molecules  with  great  relaxation  time  (for 
example,  N2); 

2)  increasing  population  of  high  levels  of  working  molecules 
through  anharmoniclty  (see  equation  (2)).  The  considerable  decline 
in  density  in  the  case  of  a  flow  in  a  nozzle  brings  about  a 
corresponding  decrease  in  the  delay  of  the  chemical  reaction. 

Under  conditions  of  a  break  between  the  vibrational  temperature 
and  the  translational  the  population  of  the  rather  high  vibrational 
levels  when  distribution  (2)  can  be  made  may  increase  when  the 
translational  temperature  decreases.  As  a  result  one  might  expect 
a  corresponding  increase  in  the  velocity  of  certain  dissociation 


m 

m 


and  exchange  reaction  as  the  temperature  declines. 

We  have  selected  as  our  subject  of  study  the  decay  of  N20  and 
the  exchange  reaction 

n2o+coh-n2+co2. 

The  following  factors  are  responsible  for  the  selection  of  these 
systems : 

1)  The  main  contributor  to  the  relaxation  rate  should  be 
level  n,  for  which  V-T-relaxatlon  is  not  yet  substantial  (i.e., 
n<n*,  where  n*  is  the  level  at  which  the  rates  of  ''-V  and  V-T- 
processes  are  compared).  Thus,  of  the  dissociation  processes  the 
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most  similar  appeared  to  be  nonadiabatic  reactions  with  predissocia¬ 
tion. 

2)  The  activation  energy  of  the  reaction  NgO+CO  according 
to  the  data  of  [7]  equals  23  Cal/mole,  i.e..  for  vibrational 

levels  with  such  energy  n<n*. 

3)  In  molecules  with  more  than  3  atoms  it  is  difficult  to 
separate  the  vibration  temperature  from  the  translational  because 
of  the  great  rate  of  the  intramolecular  V-V  exchange. 

4)  Both  of  the  selected  reactions  are  exothermal,  and  we 
can  thus  expect  an  acceleration  of  the  reaction  due  to  the  vibra¬ 
tional  energy  which  is  released. 

The  purpose  of  the  given  study  was  to  establish  the  possibility 
of  experimentally  observing  decomposition  and  exchange  reactions 
under  flow  conditions  in  nozzles  used  in  connection  with  shock 
tubes. 

Let  us  study  a  flow  of  a  mixture  of  555  N20+.95£  N2  through  a 
plane  nozzle  (height,, d  in  critical  section  1  mm,  length  20  cm, 
divergence  angle  u=12°).  To  find  the  main  characteristics  of  the 
flow  it  Is  sufficient  to  limit  ourselves  to  approximate  analytical 
relationships 

T  =  T.  (l  +  .I/’)"'  ,,  =  p,  (l  +  AT)' <  3 ) 

which  express  gas  temperature  and  density  In  terms  of  the  Mach 
number  (M),  the  ratio  of  heat  capacity,  and  the  initial  parameters 
Tq  and  pQ  at  the  nozzle  input.  In  mixtures  with  a  low  concentration 
of  a  reacting  additive  this  approximation  is  quite  satisfactory. 

The  dependence  of  the  Mach  number  along  the  nozzle  is  determined  from 
the  solution  to  the  equation 


v* 


[■hy-Mih-i'^ ^H’"' 


where  (F/F*)  is  the  ratio  of  the  current  area  of  the  cross  section 
to  the  area  in  the  critical  section. 

Let  us  write  the  equation  for  the  rate  of  change  in  the  number 
of  vibrational  quanta  falling  within  the  antisymmetrical  vibration 
of  ^OCeg)  and  the  vibration  of  nitrogen  (eA)  [4,  8] 


~  “TT^  -  A (ffl  -r  Doxp l~T~)  ~  **(*a  t1)]t 

+(fc-f»)-)fr('^r)c 

-  ■^aTra  +  ft. (t.  -r  -  t»(e.  +  1)1  -r 

^191  ~ 


Here  e?  and  e2  are  the  equilibrium  values  of  the  numbers 

HO  *  yrm  YT 

vibrational  quanta  per  one  molecule;  ta  and  Tg  are  the  times  of 

vibrational-translational  relaxation  of  molecules  A  and  B;  NA  and 

N0  -  the  densities  of  the  molecules;  6.  and  0n  -  characteristic 
b  A  0 1  B 

temperatures;  D  -  dissociation  energy;  01Q  -  probability  of 

vibrational  quanta  exchange  between  molecules  A  and  B;  ec  -  number 

of  vibrational  quanta  per  one  molecule  A  arising,  in  the  exothermal 

reaction. 


If  we  use  the  model  representing  rapid  vibrational  quanta 
exchange  [1,  8],  then  we  get  an  equation  for  eg  under  the  condition 
eg<l  taking  into  account  variable  gas  temperature,  the  chemical 
activation  A  molecule,  and  the  dissociation  of  B  molecules 

=  -  M- 'A  (~V*)]  - 

.  «-l  dTI  .  tc-tptxp(-(QA-%m  .  .  .  1  M-VbN 
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*^-*37.  a=  i-f^xe?p(--%~) 


The  relaxation  time  of  the  mixture  of  oscillators  A  and  B  is 
a  function  of  composition,  pressure,  and  the  probabilities  of 
transfer,  which  in  turn  depend  on  the  molecular  parameters  and 
temperature.  The  probability  of  V-V  and  V-T  transitions  are 
calculated  according  to  the  SSH  theory  of  [9]:  in  deriving  the 
formula  for  relaxation  time  [10]  the  condition  of  quasi-stationarity 
was  assumed  for  all  low-frequency  NgO  oscillators.  In  calculating 
the  flow  in  the  nozzle  the  temperatures  of  symmetrical  and  deformed 
vibrations  were  assumed  to  be  equal  to  the  gas  temperature. 

Under  conditions  of  an  intense  vibrational  imbalance  the 
following  dissociation  mechanism  is  used.  Vibration  excited  N2 
molecules  effectively  pump  the  antisymmetricai  type  of  N20 
vibrations,  and  the  decomposition  of  a  N20  molecules  occurs  by 
predissociation,  i.e.,  during- the  transition  of  the  molecule  into 
the  continuous  spectrum  from  the  vibrational  level  which  is  above 
the  minimal  dissociation  energy  (for  N20  Dmin=58  Cal/mole  [7]).  In 
the  simplest  case  the  dissociation  rate  (W)  under  nonequilibrium 
conditions  can  be  represented  by  the  equation  [11] 

W  -  Z;\'i,P W exp (- -5s*.) exp  {A08 -gj)}  (i  +  /,)-*  < ? ) 

Here  Z  is  the  number  of  collisions,  T^  -  the  vibrational 

temperature  of  the  antisymmetricai  N20  vibration  on  the  first  level, 

8  -  the  factor  which  considers  the  participation  of  all  vibrational 

levels  in  the  probability  of  transition  i+l->-i;  Q’  -  probability 

of  resonance  exchange  within  limits  of  antisymmetricai  vibration, 

parameter  T  determines  deviations  from  quasi-equilibrium  at 
K  _ _ 

dissociation  boundary;  Yq=0.32cx  /p/TE^;  a  -  characteristic 

—nr* 

parameter  in  exponential  potential  e  ;  y  -  reduced  mass  of 
collision;  Dfi  -  dissociation  limit  for  B  vibration  (antisymmetricai 
vibration);  ?  -  probability  of  nonadiabatic  transition.  The 
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number  of  the  level  k  is  determined  by  the  relationship 
Ek  «  ft  [El  -  (k  -  l)i£)  <  Dmin  <  EkH  ~  (ft  -f  l,  [£*,  -JfcATi 

The  decay  kinetics  of  NgO  in  their  simplest  form  are  represented 
by  the  following  scheme: 

NjO+.V-NY+O 

O  -f  NjO  — •  N*  Oj 

We  will  examine  the  effect  of  the  secondary  reaction  in  the 
two  limiting  cases: 

1)  The  rate  of  the  reaction  (2)  is  sufficiently  great,  and 
quasi-stationarity  is  maintained  with  respect  to  the  concentration 
of  0  atoms,  i.e. , 

d  IN,0 \/it  *  -2*,  INjO]  (Ml 

2)  Cooling  "freezes”  the  reaction  of  (2)  and  its  rate  is 
negligible  small.  The  first  case  can  be  realized  only  if  almost 
all  of  the  activation  energy  of  reaction  (2)  is  overcome  by  the 
vibrational  energy  N20. 

A  numerical  solution  on  ^ne  electronic  computer  was  found 
jointly  for  equation  (6),  the  equation  of  energy,  and  the  kinetic 
equation.  Under  the  condition  of  quasi-stationarity  with  respect 
to  the  reaction  (2)  we  assume  e  =D/2  (this  value  can  be  obtained 

V 

on  the  basis  of  the  results  of  [10]).  The  results  of  calculations 
for  this  case  are  shown  in  Pigs,  la  and  lb  (Pig.  la  -  mixture  of 
95J5  N2+5£  NgO,  T=1500oK,  P=10  atm,  x1/2  “  4ime  of  half-conversion 
(1);  t  -  relaxation  time  of  antisymmetrieal  N50  vibrations  under. 
V-V-equllibrium  with  N2  vibrations  (2);  x^,=T(dT/dt)  -  character¬ 
istics  cooling  time  (3);  Tg  -  vibrational  temperature  of  anti- 
symmetrical  N20  vibrations  (4)j  T  -  gas  temperature  (5);  Pig*  lb  - 
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Figure  la,  b 


We  see  that  ,'fter  x>5  cm  the  characteristic  half-conversion 

VT 

time  t1/2  becomes  less  than  t  •,  while  the  value  of  x1/2  is  such 
that  we  already  observe  a  noticeable  decomposition  in  N20  within 
the  limits  of  the  expanding  part  of  *-he  nozzle.  The  distributions 


for  the  densities  of  N2  (1),  N20  (1),  and  N20  (2)  and  the  percent 
concentration  of  N20  (3)  are  shown  in  Fig.  2;  T=1500°K,  and  PQ= 
=10  atm.  Moreover,  as  the  flow  continues  in  a  channel  with  a 


constant  cross  section  the  rate  of 


decay  exceeds  the  relaxation  rate, 
and  the  N20  molecules  continue  to 
decompose.  The  resulting  degree 
of  expansion  will  depend  on  the 
relationship  of  the  rate  of  decay 
and  translational  heating  caused 
by  the  chemical  reaction  and  V-T 
relaxation. 


Figure  3  shows  calculations 
which  do  not  consider  the  secondary 
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reactions  for  Tq*1500°K,  Pq=10  atm,  and  the  regaining  symbols  as 
in  Fig.  1.  Qualitatively  the  nature  of  the  process  is  the  same 
as  the  preceding  case,  and  the  quantitatively  differences  belong 
primarily  to  the  rate  of  the  drop  in  the  Tg  and  the  form  of  the 
Tip  (the  nonmonotonic  nature  of  x,p  on  Fig.  1  is  caused  by  the  energy 
release) . 


Interesting  are  the 
tions  Pq=1  atm,  TqOujCT; 


results  r calculating  for  initial  condi- 
(Fig.  lb)  the  line  below  the  thermal  explo¬ 
sion  limit  for  N20.  The  calculation  shows 
that  expansion  of  the  mixture  in  the 
nozzle  leads  to  an  expansion  of  N20  even 
wnen  the  theoretically  new  transition  of 
the  critical  explosion  limit  is  realized. 

Let  us  look  at  the  reaction 

N,0  -4-  CO  X,  +  CO, 

We  will  assume  that  the  reaction 
rate  is  determined  by  the  molecule 
population  on  the  vibrational  levels  whose 
sum  energy  exceeds  activation  energy  E= 

=23  Cal/mole  [7] 


if-zpS.v.S.v* 

fmo  m 

where  N  and  N  are  the  populations  on  levels  m  and  n,  determined 
n  m 

from  che  Trinor  distribution  (see  (20).  Here  En+Em>E,  P=10  :  is 
the  steric  factor.  The  calculations  show  that  when  tho  mixture  of 
N20+C0  expands  in  a  nozzle  conditions  may  be  created  under  which 
tc<tVT*  i.e.,  an  exchange  reaction  may  occur  during  tim?s  smaller 
than  the  relaxation  t.  me .  If  the  release  of  energy  occurs  primarily 
in  vibrational  degrees  of  freedom  of  N2  and  CO,  then  this  will  lead 
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to  a  further  acceleration  in  the  reaction.  Also  advantageous  from 
the  standpoint  of  maintaining  equi]ibrlum  is  the  fact  that  the 
reaction  products  have,  generally  speaking,  greater  relaxation 
times  than  the  original  reagents. 


Despite  the  approximate  nature  of  the  calculation  (resulting 
primarily  from  the  fact  that  there  does  not  exist  a  sufficiently 
reliable  dissociation  theory  under  nonequilibrium  conditions),  we 
believe  that  the  experimental  observation  of  the  decomposition 
during  cooling  in  a  supersonic  expanding  flow  is  reliable.  Moreover 
the  possibility  of  exchange  reactions  (N20+C0,  KI+F2,  and  others) 
occurring  during  times  less  .than  the  V-T  relaxation  times  is 
extremely  interesting. 
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